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Abstract: We present a new method for solving the BFKL evolution applicable at both 
leading and next-to-leading logarithmic accuracy, and tailored to the study of QCD multi- 
jet events at colliders. We utilise this to discuss corrections to the standard analysis. There 
are known, large corrections from energy and momentum conservation. We show that, 
despite claims to the contrary in the literature, these are unrelated to the next-to-leading 
logarithmic corrections to the evolution kernel. 
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1. Introduction 

One of the many immediate challenges for QCD is to provide a reliable description of the 
multiple hard jet environment of the LHC. Besides posing a very interesting problem in 
itself, the QCD dynamics will provide signals similar to that of many sources of physics 
beyond the standard model, and so is very important to understand in detail. An intriguing 
alternative to the standard approach of calculating the production rate of a few hard 
partons by fixed order perturbation theory is to use the framework arising from the BFKL 
(Balitskii-Fadin-Kuraev-Lipatov) equation to calculate the emission of gluons (and quarks 
at next-leading logarithmic accuracy) from the evolution of an effective, Reggeized gluon 
(Reggeon) propagator. The starting point here is the observation that for e.g. 2 —* 2, 2 — > 
3, . . . gluon scattering, the partonic cross section in the limit where the rapidity span 
of the two leading gluons is increased is dominated by the contribution from Feynman 
diagrams with a t-channel gluon exchange. This t-channel gluon is then evolved according 
to the BFKL equation, and will emit partons accordingly. Starting from the 2 — > 2- 
gluon exchange, the 2 — > 2 + n gluon scattering process can be calculated in the limit of 
large rapidity spans A, thanks to the Regge factorisation of the colour octet exchange. 
Obviously, this means that the formalism is relevant only if there is sufficient energy at 
colliders to have multiple emissions spanning large rapidity intervals. In this high energy 
limit, the partonic cross section (p' a ,p' b — > p a , {Pi},Pb) factorises as 

d°(Pa,Pb) = Ta(Pa) /(Pa, ~P&, A) r 6 (p 6 ), (1.1) 

where p a ,Pb is the momentum of the partons furthest apart in rapidity, T a ^ are the process 
dependent impact factors, describing here the gluon-gluon-Reggeised-gluon-coupling, and 
f(Pa, — P&, A) is the process independent gluon Green's function, which is the object evolv- 
ing according to the BFKL equation. We use bold face vectors to denote the transverse 
components. This is the same gluon Green's function that would enter calculations of the 
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small- x evolution of parton density functions. The evolution of the gluon Green's function 
can of course be discussed generally without reference to specific impact factors, and so the 
conclusions of this study should impact all applications of the BFKL formalism for colour 
octet or inelastic studies. In Eq. ( |1.1| ), the partonic cross section has been integrated over 
any number of gluons emitted in this evolution, so only the dependence on the momenta of 
the leading jets have been retained. The BFKL approach should give a good approximation 
to the full calculation, when the hard jets are well separated in rapidity. 

The BFKL equation governing the evolution of the t-channel Reggeised gluon was 
solved iteratively in Ref. [1-3] using the evolution at leading logarithmic accuracy [4-6], 
and more recently in Ref. [7,8] at full next-to-leading logarithmic accuracy [9,10]. The 
present paper has three purposes. Firstly, to present a new method for obtaining the QCD 
evolution according to the BFKL equation. Secondly, to discuss the implications of this 
new formalism on our understanding of the sources of corrections. And thirdly, to announce 
the availability of a computer code that calculates the exclusive multi-parton production 
rate expected at the LHC or Tevatron according to the BFKL evolution. 



2. Iterative Solution of the BFKL Equation 

Our starting point for the discussion is the fully inclusive BFKL equation describing the 
evolution of the gluon Green's function /(k a ,kf,, A) 



co UKM) = 5 (2+2e) (k a - k 6 ) + j d 2+2s k £ e (k a , k + k a ) / w (k + k a , k 6 ) , (2.1) 

where w is the Mellin-conjugated variable to A, and the BFKL kernel /C e (kj, kj) is presently 
known to next-to-leading logarithmic accuracy, where the logarithm is ln(sy/|kj||kj|), and 
Sjj the invariant mass of partons i and j. The solution to the evolution from a momentum 
kb at rapidity y^ to k a at rapidity y a according to this integral equation can be written on 
the form 

oo „ 

/(k a ,k 6 ,A) = ^ / dV n Fn, 
n=0 

JdV n =(jf[ dki jHdj/i fdy 2 ■ ■ ■ f V Jdy n ^ tf 2 > (k a + £ k, - k^ (2.2) 



n 



Tn = IJe^dfc-i-*) y( qi ,q J+ i) e 



,w(q n +i)(l/ n -l/n+l) 



\i=l 



with yo = y a = A, y^ = y n +i = 0, = k a + ^2]Z\ k^, and, crucially, the real production 
vertices V(qj,qj+i) and trajectories u>(qi)(yi-i — yi) are regularised and finite at LL [1-3] 
and NLL [7, 8] in order to facilitate a direct numerical evaluation. The correctness of the 
procedure at NLL was proved in Ref. [11] by comparing to analytic results for the evolution 
in Af =4 Super Yang-Mills theory, since the analytic solution to the BFKL equation is only 
known for conformal invariant theories (please see Ref. [12] for a discussion of the analytic 
methods at NLL). 
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Figure 1: The 2 — > 2 + n gluon scattering process described using Regge factorisation, with the 
initial state at the bottom. The shaded blobs are the gluon-gluon-Reggeon impact factors T a ,b, 
and the hatched blobs are the (regularised) gluon-Reggeon-Reggeon vertices. The Reggeized gluon 
(Reggeon) propagators are marked with zigzag lines. Gluon emission is generated in the rapidity 
span between the impact factors by the evolution described by the BFKL equation of the Reggeized 
gluon. At NLL the vertices can emit one or two gluons, or a quark-anti-quark pair. 



The formulation of the solution to the BFKL evolution in Eq. (|2.2| ) allows in principle 
for the calculation of the fully exclusive 2— >2 + n cross section, i.e. the expression of the 
differential cross section in Eq. (|1.1| ) as da(p a , {pi},Pb)- It is reassuring that this picture 
is fully consistent with recent results on the multi-Regge form of QCD amplitudes at the 
next-to-leading logarithmic level [13-15]. 

3. Direct Solution of the BFKL Evolution 

Starting from Eq. ( |2,2j ) it is possible to construct a new method of evolving according to 
the BFKL kernel by first performing a simple change of variables in the nested integration 
over rapidities, followed by the introduction of another integral and delta functional, in 
order to make the integration over rapidity separations independent 

<-yo rdyi rdy n 



J d Vl J dy 2 ■■■ J dy n f ]J ^Xvi-i-ih) J e ^+i)(y n -yn+i) (3 -Q 

„A /-A-j/n r^~yn 2/2 ( n \ 

J^d5y n J^ d5y n ^---^ dS Vl I JJ^(«w)«wJ (3.2) 

roo roo roo "+1 "+1 

/ d5y n+1 / dSy n --- dS Vl «5(A - V^) TT (3.3) 
Jo Jo Jo i=1 i=1 



so the gluon Green's function /(k a , k^, A) in Eq. ( |2.2| ) can be written as 

00 „ 

/(k a ,k 6 ,A)= J] / dV n ^ n , 
n=0 

/ n roo \ roo / n \ / n+1 \ 

dV n =\JlJ dk * J o d6 y*J J Q d6 Vn+l ^ ( ka + ^ ^ " ^ J 5 [ A ~ E S V* J 



(3.4) 



Tl 



Tn= m e ^)^ V(qi,q i+1 ) e> 



^(qn + l)<5j/n 



+ 1 



\i=l 
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Using this, we see that /(k a , k&, A) is simply the value at A = Y27=i ^Vi °^ ^ ne P^ase space 
integral of the product of real-emission vertices V(qj,qj+i) at rapidity yi = Y^j=i^Vj 
connected with Regge factors e UJ ^ i ^ Syi describing the probability of no (resolved) emission 
between two adjacent (in rapidity) vertices 1 . In this case, it is clear that the virtual and 
unresolved corrections encoded in cj(qj) lead to a decrease with increasing A, while any 
increase in /(k a ,k&, A) is due to the integration over phase space of the resolved emission 
from the vertices V(qi,qi+i). The BFKL evolution can then be found by the following 
algorithm :- 

1. Choose a random number of vertices for the evolution, n > 

2. Generate a set {kj}j = i n of transverse momenta (the outgoing momenta are {— kj}j = i r .. jn ) 

3. Calculate the corresponding set of trajectories {w(qi)}i=i,..., n +i, and vertex factors 
{V{q.i,q.i+i)}i=i,...,n, q« = K + Ya=i^i 

4. Generate the inter- vertex rapidity separations {dyi} according to the distributions 

e uj(m)Syi 

5. Calculate the corresponding A = Y17=l ^Vi anc ^ return Y\7=i Vi^ki 

6. Repeat until required Monte Carlo accuracy is obtained 

This algorithm is vastly superior to the immediate Monte Carlo implementation of the 
phase space integrals in Eq. (|2.2j). This is because that instead of calculating, for a given 
A, the contribution from any number of emissions and any momentum configuration of 
these, this new method calculates a representative rapidity span for a given number of 
emissions and their configurations of momenta. Furthermore, this method is explicitly 
symmetric in k a <-> kj, and evolution direction (increasing or decreasing rapidities), and it 
offers a more direct way of implementing the necessary study at colliders of any length of 
the rapidity interval A, since this is automatically achieved by just a single sum over any 
number of emissions. Unweighting of the Monte Carlo is also significantly more efficient, 
since part of the integrand in the previous formulation contributing to the variance of the 
integrand is now used to determine a representative value for A. We stress that both 
approaches provide the same correct solution to the BFKL evolution. 

4. Phase space restrictions from energy and momentum conservation 

Besides leading to a more efficient implementation of the BFKL evolution, the new for- 
mulation of the solution to the BFKL evolution is also useful for discussing sources of 
corrections to the standard BFKL analysis. The phase space restrictions from energy and 
momentum conservation can be implemented simply by adding a step in the algorithm 
rejecting events prohibited by phase space considerations. We see that the discussion of 
the region of integration for real emission is completely independent of the logarithmic 

1 This interpretation holds for w(qi) < 0, which can always be fulfilled in the regularisation procedure of 
Ref. [1-3] at LL and Ref. [7,8] at NLL. 
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accuracy of the evolution kernel expressed in terms of the vertices V(c[i, qi+i) and Regge 
trajectories w(qj). This clearly shows that, despite what is often claimed in the literature 
(see e.g. Ref. [16] and references therein), phase space restrictions in terms of e.g. energy 
and momentum conservation are completely independent of the logarithmic accuracy to 
which the BFKL evolution is performed. Specifically, the NLL corrections to the evolution 
kernel do not implement energy and momentum conservation. Furthermore, the standard 
solution of the LL and guestimate of the NLL evolution based on a Mellin transform in 
the transverse momentum of the kernel will always fail to take such considerations into 
account, by its very nature as fully inclusive both in number of emissions and as an inte- 
gral to infinity of the transverse momentum. This is simply a result of the fact that the 
BFKL evolution is local, while energy and momentum conservation depends on the full 
final state configuration. Such considerations should be implemented by modifications to 
the evolution beyond the discussion of the logarithmic accuracy of the evolution. Please 
note that this discussion of energy and momentum conservation goes beyond the discussion 
of longitudinal momentum in collinear splittings at small- x , implemented by ensuring the 
vanishing of the first moment (see e.g. Ref. [17]). 

While it is not our job to guess the cause of the confusion in the literature over 
the role of energy and momentum conservation and higher logarithmic corrections to the 
evolution kernel, it is perhaps beneficial to discuss the differences and similarities between 
the multi-particle generating colour-octet exchange and the colour-singlet exchange relevant 
for diffractive studies. The evolution in rapidity A of the gluon Green's function for 
both cases is described by a BFKL equation. In the case of diffractive 2^2 processes, 
the rapidity span of the BFKL evolution is given by A ~ lns/so, when the Regge scale 
so = | k a 1 1 kft j is the product of the transverse momentum of the two (massless) scattered 
particles. However, the impact factors depend on so only at NLL accuracy, and so it can 
be argued that only at NLL accuracy does the prediction gain a correct dependence on the 
centre of mass energy (although of course the evolution can (and should!) be discussed 
without reference to specific impact factors). While this is true, care has to be taken when 
discussing instead colour-octet exchange with multiple emissions. First of all, there is 
no one-to-one correspondence between the rapidity span of the evolution and the centre of 
mass energy, although obviously the centre of mass energy tends to increase with increasing 
rapidity span. Note specifically that if one sets s = soe A with the Regge scale sq = |k a ||k&|, 
then any additional emission from the BFKL evolution is kinematically excluded. The 
error 2 in ignoring the contribution from the BFKL emission to the centre of mass energy 
is significant [18,19]. Secondly, any constraint on the real emission will significantly lower 
the expected BFKL signatures, e.g. the expected rise in Fi at small-x, and the increasing 
jet azimuthal decorrelation with rapidity at hadron colliders. In fact, the effect of correctly 
implementing energy and momentum conservation on top of the LL BFKL evolution has a 
larger impact on the jet azimuthal decorrelation even at LHC energies than the inclusion 
of the NLL corrections to the evolution kernel (compare e.g. Fig. 10 of Ref. [8] with Fig. 5 
of Ref. [3]). 

2 It is often argued that this correction is logarithmically subleading. However, as we have demonstrated 
this is not part of the subleading corrections to the evolution kernel. 
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5. Conclusions 



We have presented a new method for solving the BFKL evolution very effectively for the 
study of multi-partonic final states at hadron colliders. We have furthermore demonstrated 
that the discussion of energy and momentum conservation and phase space constraints 
of the evolution in general is completely separate to the discussion of the logarithmic 
accuracy of the evolution kernel, contrary to the claims found in the literature. We have 
demonstrated how energy and momentum conservation can be implemented. 

A program implementing the new algorithm for the BFKL evolution in the case of 
multi-jet production at Tevatron and LHC energies is available at the URL 
http : //www . hep . phy . cam . ac . uk/~andersen/BFKL 

The current version implements energy and momentum conservation, and the BFKL evo- 
lution kernel at LL supplemented by the running coupling terms from NLL. 
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